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Abstract 

A model for charge transport in undoped, photo-excited semiconductor superlattices, which 
includes the dependence of the electron-hole recombination on the electric field and on the photo- 
excitation intensity through the field-dependent recombination coefficient, is proposed and ana- 
lyzed. Under dc voltage bias and high photo-excitation intensities, there appear self-sustained 
oscillations of the current due to a repeated homogeneous nucleation of a number of charge dipole 
waves inside the superlattice. In contrast to the case of a constant recombination coefficient, nu- 
cleated dipole waves can split for a field-dependent recombination coefficient in two oppositely 
moving dipoles. The key for understanding these unusual properties is that these superlattices 
have a unique static electric-field domain. At the same time, their dynamical behavior is akin to 
the one of an extended excitable system: an appropriate finite disturbance of the unique stable 
fixed point may cause a large excursion in phase space before returning to the stable state and 
trigger pulses and wave trains. The voltage bias constraint causes new waves to be nucleated when 
old ones reach the contact. 

PACS numbers: 73.63.Hs, 05.45.-a, 73.50.Fq, 72.20.Ht 
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I. INTRODUCTION 



Nonlinear charge transport in weakly coupled, undoped, photo-excited, type-I semicon- 
ductor superlattices (SLs) is well described by spatially discrete drift-diffusion equations.-^ 
As in the much better known case of doped SLs, nonlinear phenomena include formation 
and dynamics of electric-field domains, self- sustained oscillations of the current through 
voltage-biased SLs, chaos, etc. Experimentally, the formation of static electric-field domains 
in undoped, photo-excited SLs was already reported many years ago. 3 The first experimental 
observation of dynamical aspects of domain formation in undoped, photo-excited SLs such 
as self-sustained oscillations of the photo-current were reported by Kwok et alA Due to the 
excitation condition, the oscillations were damped. Subsequently, undamped self-sustained 
oscillations of the photo-current in undoped SLs were observed for a type-II GaAs/AlAs 5 
and for a direct-gap GaAs/AlAs SL.- Tomlinson et all- reported the detection of undamped 
photo-current oscillations in an undoped GaAs/Al .3Gao.7As SL, where the transport is gov- 
erned by resonant tunnelling between V states. The evolution from a static state at low 
carrier densities to an oscillating state at higher carrier densities was demonstrated in an 
undoped, photo-excited SL by increasing the photo-excitation intensity. 8 An investigation of 
the bifurcation diagrams for undoped, photo-excited SLs showed the existence of a transition 
between periodic and chaotic oscillations. 9 For a detailed review of the nonlinear static and 
dynamical properties of doped and undoped superlattices, see Ref . [10 . 



Previous theoretical studies of undoped photo-excited SLs, including studies of bifurcation 
and phase diagrams,™ are based on a discrete drift model having a constant recombination 
coefficient.- So far, there have been no reports on considering field-dependent recombination 
or the fact that the time scale of the electron-hole dynamics depends strongly on the optical 
excitation intensity. However, the consequences of including these effects for the dynamics 
of electric-field domains can be striking. In this work, we incorporate into the previously 
studied discrete model the dependence of the electron-hole recombination on the electric field 
and on the photo-excitation intensity using a straightforward model that takes into account 
the overlap integral between the electron and hole wave functions. The field-dependent 
recombination coefficient decreases with increasing electric field, which has far reaching 
consequences. 

At high photo-excitation intensities, it is possible to find only one stable electric-field 
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domain, not two as in the case of a constant recombination coefficient.- In this case, self- 
sustained oscillations of the current (SSOC) may appear under dc voltage bias. The field 
profile during SSOC can exhibit nucleation of dipole waves inside the sample, the splitting of 
one wave into two, and the motion of the resulting waves in opposite directions. These dipole 
waves resemble the pulses in excitable reaction-diffusion systems such as the FitzHugh- 
Nagumo model for nerve conduct ion^^^^ and are quite different from field profiles for 
a constant recombination coefficient.-^^ In an excitable dynamical system, an appropriate 
finite disturbance of the unique stable fixed point may cause a large excursion in phase 
space before returning to the stable state. When diffusion is added, the resulting reaction- 
diffusion system may support wave fronts, pulses, and wave trains. It is also possible to 
find SL configurations at high photo-excitation intensities for which there exist no stable 
electric-field domains. In these cases, there are SSOC, whose corresponding field profiles are 
wave trains, comprising a periodic succession of dipole waves. These cases are similar to 
wave trains in oscillatory media such as those appearing in the FitzHugh-Nagumo model in 
the presence of a sufficiently large external current.— 1 ^ 

II. MODEL EQUATIONS 

The equations governing nonlinear charge transport in weakly coupled, undoped, photo- 
excited, type-I SLs are 

e(Fi-Fi_i) =e{m-pi) , i = l,...,N, (1) 
dF- 

£ ^ + ^ +1 = J(t) ' 2 = ' 1 '---> iV ' ( 2 ) 
d ^ = 1 (I)-r{F l J)n lPl , i = l,...,tf. (3) 

In these equations, the tunneling current densities between the quantum wells (QWs) as 
well as between the SL and the contact regions are 



Ji-H+i = -j ~ eDi(Fi) l+1 p , (4) 



Jo^i — v Fo , Jn^n+i — "tt ® Fn ■ (5) 



The voltage bias condition is 



1 N V 

i=0 



Here — e < 0, e, a, Np, —Fi, rii and pi denote the electron charge, the average permittivity, 
the conductivity of the injecting contact, the doping density of the collecting contact, the 
average electric field, as well as the two-dimensional electron and hole densities of the ith 
period of the SL, respectively. Equation ([T]) corresponds to the averaged Poisson equation. 
Equation (j2J) denotes Ampere's law: the total current density J(t) equals the sum of the 
displacement current density and J^+i, the electron tunneling current density across the 
iih barrier that separates the quantum wells i and i + 1. Charge continuity is obtained by 
differentiating Eq. ([1]) with respect to time and using Eq. ([2]) in the result. Tunneling of 
holes is neglected so that only photogeneration and recombination of holes with electrons 
enter into Eq. For high temperatures, i.e. ksT ^> nh 2 [n]/m*, where ks denotes 

Boltzmann's constant, T the temperature, m* the effective electron mass, and [n] the order 
of magnitude of the electron density, the tunneling current is given by Eq. (j3J), in which v(Fi) 



10 and / = L w + L B is the length 



and D(Fi) are functions of the electric field given in Ref. 
of one SL period (L w and Lb denote the individual widths of the quantum well and barrier, 
respectively). Even for lower temperatures, the qualitative behavior of the solutions of the 
discrete drift-diffusion model is similar to one of the more general tunneling current models 



described in Ref. 



1CJ. For D = and a constant recombination coefficient r, Eqs. (JTJ) — (| 



describe the well known discrete drift model introduced in Ref. The indices i — and 
i = N + 1 represent the SL injecting and collecting contacts, and Eq. (j2J) holds for them 
with the phenomenological currents given by Eq. (jSJ). The total current density follows from 
the voltage bias condition in Eq. (jHJ) 

A' 

,v + " 



i=0 



Photogeneration and recombination are given by ■j(I) = I a^ihjj^) L w /(hu exc ) and 

1 J Wire) J exp(^) - 1 V } 

respectively. Here J, iv exc , n re f, n in w 'jn re fL w /c, and c denote the photo-excitation inten- 
sity, the frequency of the exciting photon, the refractive index, the intrinsic carrier density, 
and the speed of light, respectively. « 2 d and « 3 d correspond to the two- (2D) and three- 
dimensional (3D) absorption coefficients. The 2D absorption coefficient is proportional to 
the square of the modulus of the electron-hole overlap integral for a constant electric field 



4 



F (cf. Ref. Il6|) 



a 2D (u,F) 



of I 8{E h + E g -hw + E n e {F) + E™{F)) dE h 
o 

1/2 



(9) 



x 



^ e {z,F)^ h {z,F)d~z 



-1/2 



where E g denotes the energy of the bandgap at the T point and as well as ^ solve 
the stationary Schrodinger equation inside one SL period, [—1/2, 1/2], for the electrons and 
holes, respectively. In this equation, the electric field F is considered to be constant and 
ty n! h(±(L w /2 + l p ±)) = 0, where the penetration length l p ± solves the cubic equation 



2m* 



V-eF[l p± T 



L ,, 



E„ 



n,h 



and therefore l p ± depends self-consistently on the eigenvalue E n ^. For a fixed value of /, 
the recombination coefficient decreases with increasing electric-field strength F, as depicted 
in Fig. [H 



III. DIMENSIONLESS EQUATIONS 

To study the model described by Eqs. ([I])-©, it is convenient to render them dimen- 
sionless. For a fixed value of the photo-excitation intensity / and a constant electric field 




FIG. 1: (Color online) Recombination coefficient r(F,I)/r(0,I) vs electric field F for the four 
configurations listed in Table HI 
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Fi = F, the stationary solution of Eq. is rij = Pi — ^/ / y(I)/r(F, I). We use the maximum 
values [n] = [p] = a/7(/)/V(0, I) to define typical values for and pj. The field F M at the 
maximum of the drift velocity is a typical value of the field when there are SSOC Therefore, 
we adopt it as the field unit [F] = Fm- Similarly, [v] = vm and, therefore, [J] = e[n]%/l 
and [D] = Ivm- There are two possible time scales, the first one being tp — eF M /[J], which 
balances Maxwell's displacement current with the current density in Eq. (j2J), and the second 
one t n = \p\/l{F) = l/y / 7(/)r(0, 1), which balances both sides of Eq. ([3]). It is reasonable 
to choose the time unit as the longer of the two times tp and t n . We have chosen four rep- 
resentative GaAs/ALcGai-ajAs SL configurations with 10 nm wells and 4 nm barriers. The 
configurations are I: x — 0.25 and / = 120.5 kW/cm 2 ; II: x = 0.25 and / = 302.69 kW/cm 2 ; 
III: x = 0.3 and / = 479.735 kW/cm 2 ; and IV: x = 1.0 and / = 30.27 kW/cm 2 . We assume 
a circular cross section with a diameter of 160 /im, a photo-excitation intensity of 60 mW, 
a wavelength of 413 nm, and four different beam diameters yielding the previously listed 
values of the laser intensity. For most of the cases listed in Table [H t n > tp, and therefore 
we choose [t] = t n . All scaled variables are listed in Tabled! 

We now rewrite the model equations using dimensionless variables by defining hi = ni/[n], 
t = t/[t], where [n], [t], etc. are the scales defined above and specified in Table [0 
Omitting hats over the variables, the dimensionless system of equations corresponding to 

TABLE I: Units used to achieve a set of equations with dimensionless variables. 
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Eqs. (P)-© read 



Fi — = (m — Pi)v , (10) 
dF 

6 ~di + UiVi - Fl) ~~ D{Fl) {ni+1 ~ Hl) = J(t) ' 

d J^ = \-r{F l )n lVll (12) 

o F + 5 — = J , aan N F N + 5 — ^- — J , (13) 
at at 

1 " 

i=0 

i=0 

In these equations, there are four dimensionless parameters, 

eF M /r(0,J) _ e[n] _ [n] _ yftij) 



ev M eF M N D N Dy /r(0,I) 

and the dimensionless conductivity a, which is here simply denoted by a. The values of 
the dimensionless parameters for the four SL configurations given in Table d are listed in 
Table HT1 It is interesting to note that r(0, I) oc I~ 2 and 7(f) oc / so that v oc I 3 ^ 2 and 
5 oc J -2 , i.e. v increases with photo-excitation intensity, whereas 6 decreases. High photo- 
excitation intensities imply that v is large and 5 small, whereas the opposite holds for low 
photo-excitation intensities. 

TABLE II: Numerical values of the dimensionless parameters a, 5, v and a for the four superlattice 
configurations listed in Table HI 



config. a 5 v a 

M eF M lr(0,I) e[n] a F M l 
Nd eVM sFm evM N 

I 0.9713 0.0037 8.446 1.0834 



II 0.6944 5.91 x 1(T 4 33.6228 1.0523 

III 0.9779 4.05 x 1(T 5 70.507 1.1074 

IV 1.5754 22.3714 1.2721 1.2959 
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It is interesting to depict the phase plane corresponding to spatially uniform solutions of 
Eqs. (jlDjl-ipj) with n i =p i =p,F i = F 

5^ = J-pv(F) } f t =l-r(F)p 2 . (17) 

Generically and for a fixed value of J, the nullclines v(F)p = J and r(F)p 2 = 1 intersect 
in one or three fixed points, depending on the Al content x in the barriers. At these fixed 
points, 

3(F) = J, J'(^) = 4=- (18) 
s/r(F) 

The function j(F) is depicted in Fig. [2] for the SL configurations listed in Table [H 

The calculations for arbitrary values of Aluminum content x show that for 0.45 < x < 1, 
there are three fixed points of the system in Eq. (|T7|) . one on each of the three branches 
of p = J/v(F), and the ratio of (j max - j min ) to the average current (j max + j min )/2 is 
sufficiently large. As we shall see later, some nonlinear phenomena occurring in these SLs 
are quite similar to the ones observed in doped SLs: static electric-field domains with domain 
walls joining the stable branches of p = J/v(F), SSOC due to the recycling of pulses formed 
by two moving domain walls having a high-field region between them, etc. For < x < 0.45, 
j(F) is either increasing for positive F (for < x < 0.25) or the ratio of (j 

max Jmin) to 

the average current (j max + jmin)/2 is small (for 0.25 < x < 0.45). For < x < 0.25, 
there is a unique fixed point at F = F#, which, for an appropriate value of J, may be 




F 



FIG. 2: (Color online) Local current density j vs electric field F. For a fixed value of the total 
current density J, there may be one or three zeros of j(F)-J depending on the Al content x. 
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located on any of the three branches of p = J/v(F). If the fixed point is located on one of 
the two stable branches of p = J/v(F) for which v(F) has a positive slope, the dynamical 
system of Eq. (IP7|) is excitable, whereas it is oscillatory if the fixed point is located on the 
second branch of p = J/v(F) with v'(F # ) < —25r{F if ) < 0. In an excitable dynamical 
system, an appropriate finite disturbance of the unique stable fixed point may cause a large 
excursion in phase space before returning to the stable state. When diffusion is added, the 
resulting excitable react ion- diffusion system may support a variety of wave fronts, pulses 
and wave trains.— A dynamical system having an unstable fixed point and a stable 
limit cycle around it is oscillatory. Again in the presence of diffusion, oscillatory systems 
may support different spatio-temporal patterns.— An undoped photo-excited SL with 
excitable or oscillatory dynamics exhibits quite unusual phenomena. Under dc current bias, 
it is possible to have pulses moving to the right or to the left and periodic wave trains. 
Under dc voltage bias, these pulses and wave trains may give rise to a variety of SSCOs. 
Similar phenomena are observed in the case 0.25 < x < 0.45 for which j(F) has a shallow 
valley for a narrow interval of current densities. 

IV. DC VOLTAGE-BIASED SUPERLATTICE FOR SMALL PHOTO- 
EXCITATION INTENSITIES 

The behavior of a dc voltage-biased SL is quite different depending on its Al content and 
photo-excitation intensity. For an Al content smaller than 45%, j(F)-J in Eq. ffT8]) has a 
single zero for any value of J, and the only stable states of the SL are stationary ones unless 
the photo-excitation intensity is sufficiently large (cf. next section). Let us assume that 
the Al content is larger than 45% so that j(F)-J in Eq. (1181) may have three zeros for an 
appropriate range of J values. In this case, the undoped SL behaves similarly to an n-doped 
SL.— The most interesting limit is that of small photo-excitation intensity, i. e., (5>1. 

For the time scale r = t/8, Eqs. ( 11 II) and (fl2|) can be rewritten as 




(19) 



(20) 
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In the limit of large photo-excitation intensities, S <C 1, Eq. (|20|) indicates that the p, do 
not depend on r. In this case, Eq. ffl9l) corresponds to a SL doped with a density and we 
may expect phenomena similar to the ones observed in an n-doped SL. In the opposite limit 
of small photo-excitation intensities, C 1, r = 5~ 1 t is a slow scale. The Pi are functions 
of Fi — Fi-i determined by solving Eq. (1201) with a zero left hand side 



Pi 



'p. _ p. 
i 1 i— 1 



r(Fi) 



F - F 1 
2^ 



(21) 
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FIG. 3: (Color online) Current vs time (top) and field distribution vs time (bottom) displaying 
SSCOs due to monopole recycling at the injecting contact. The injecting contact conductivity and 
the bias are a = 1.2959 and V = 2.0416236, respectively. The other parameter values correspond 
to configuration IV in Table UJ and are N = 99, 5 = 22.35636, a 2 = 1.3148, v = 1.27214, and Al 
content x = 1. 
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The resulting equation for F, is then obtained by inserting Eq. ([ST]) into Eq. (|T9l) 




+ 



+ 1+1 £ D(F t ) . (22) 

This equation is similar to the one describing the electric field in a doped SL, but now there 
are drift and diffusion terms which are nonlinear in the differences Fi — F^\. Under dc 
voltage bias, there are SSOC mediated by pulses of the electric field. The current density 
varies on a slow time scale, whereas the electric-field profile consists of a varying number of 
wave fronts joining the stable constant solutions of Eq. (122]) at the instantaneous value of the 
current.— Depending on the conductivity of the injecting contact a, there are different types 
of SSOC due to the periodic generation of dipoles or monopoles at the injecting contact. If 
the curve aF intersects the bulk current-field characteristic curve j(F) before its maximum 
(cf. Fig. [2]), SSCOs due to recycling and motion of charge monopole waves (moving charge 
accumulation layers) appear as shown in Fig. If oF intersects j(F) after its maximum 
(cf. Fig. |2]), SSCOs due to recycling and motion of dipole waves are obtained, as depicted 
in Fig. HI We have indicated in Figs. [3] and H] whether there are one or three uniform and 
time-independent (constant) solutions of Eq. (|17p . solving pv(F) = J and r(F)p 2 = 1 for 
the instantaneous value of the current density J = J(t) during the SSOC. 

V. DC VOLTAGE-BIASED SUPERLATTICE FOR LARGE PHOTO- 
EXCITATION INTENSITIES 

If < x < 0.45, the phase plane in Eq. (fTTj) may have only one fixed point located on any 
branch of the nullcline p = J/v(F). For small photo-excitation intensities, the only stable 
state is a stationary one. However, for sufficiently large photo-excitation intensities, an 
infinite, dc current-biased SL may exhibit pulses moving downstream or upstream and also 
wave trains moving downstream.— The counterpart of these stable solutions for a dc voltage- 
biased SL is very interesting and different from anything observed in an n-doped SL. Our 
simulations correspond to SLs with different Al contents, whose current-field characteristics 
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FIG. 4: (Color online) Current vs time (top) and field distribution vs time (bottom) displaying 
SSCOs due to dipole recycling at the injecting contact. The injecting contact conductivity and the 
bias are a = 0.5625 and V = 2.05148375, respectively. The other parameter values are the same 
as the ones used for Fig. [3] corresponding to configuration IV in Table HI 

j(F) and injecting contact curve j = oF are shown in Fig. [2j In all cases, SSCOs appear 
for an average bias roughly in the region of negative differential resistance (NDR), where 
j'(F) < (e.g., (T=0.5625 in Fig. [2D- 

When the conductivity of the injecting contact is such that <jF intersects j(F) near 
the maximum thereof, it is possible to have SSCOs that are quite different from the ones 
appearing in n-doped SLs. For an injecting contact conductivity a = 1.05231 (cf. Fig. [2J), 
pulses (charge dipoles) may be triggered at the injecting contact, move toward the receiving 
contact, and cause SSCOs as shown in Fig. [51 For all the instantaneous values of J(t) 
during these SSCOs, there is only one constant solution of Eq. ( TT7I) : most of the times 
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FIG. 5: (Color online) Current vs time (top) and field distribution vs time (bottom) displaying 
dipole-mediated bulk SSCOs for an injecting contact conductivity a = 1.05231 and x = 0.25 
(cf. Fig. [2|). These oscillations correspond to having a finite wave train. The other parameter 
values correspond to configuration II in Table U and are N = 99, v = 33.62275, 5 = 5.91 x 10~ 4 , 
V = 1.22732, and a 2 = 0.73075. 



this solution is on the second NDR branch of j{F). Only when J(t) is near its maximum 
value, the constant solution is on the third branch of j(F). For a constant current density, 
the constant solution of Eq. (ITT)) is on the NDR branch, which implies that the system is 
oscillatory and that periodic wave trains are possible. The realization of wave trains for a 
long dc-biased SL are displayed in Fig. [6j at any time during SSOC, there are only two 
fully developed pulses present in this SL. These pulses experience variations in their shape 
and velocity when they are generated or arrive at the contacts, but longer SLs allow for 
realizations of wave trains, in which more pulses exist simultaneously inside the SL. 

In the previous example, pulses always move downstream, from left to right. For slightly 
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FIG. 6: (Color online) (a) Numerically obtained field profile of a pulse moving with positive speed 
followed by a wave train when there is only one critical point in the phase plane, (b) Phase plane 
showing the nullclines and the motion of the 350th QW as the pulse traverses it. The total current 
density is J = 1.0165 and the other parameters correspond to configuration I in Table HI 



larger conductivity of the injecting contact (cr = 1.107 in Fig. |2J), Fig. [7| shows that two 
pulses are formed inside the SL and move with opposite velocities toward the contacts. 
These pulses are similar to the ones constructed above for the case of dc current bias (with 
positive or negative velocity), except that the current changes slowly with time during the 
self-oscillation and the pulses accommodate their form to the instantaneous value of the 
current. Note that the pulses are triggered inside the SL, not at the injecting contact. A 
pulse moving with positive speed has a long trailing region, in which the field increases as 
we move away from the pulse. However, there is a depletion layer near the injecting contact, 
in which the field decreases as we move away from the injecting contact. In a long, but finite 
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FIG. 7: (Color online) Current vs time (top) and field distribution vs time (bottom) displaying 
dipole-mediated bulk SSCOs for an injecting contact conductivity a = 1.107 (cf. Fig. [2]). The 
other parameter values correspond to configuration III in Table Hand are N = 61, v = 70.50677, 
5 = 4.069 x 10~ 4 , V = 1.269848, and a 2 = 1.0829727. 

SL, a local maximum is formed inside the SL, when the decreasing field near the contact 
meets the increasing field in the trailing region of the exiting pulse with positive speed. The 
current increases as the exiting pulse is absorbed by the receiving contact, until it surpasses 
a critical value. In this case, the local maximum of the field profile inside the SL is split, and 
two new pulses are created. The pulse closer to the injecting contact moves toward it with 
negative speed whereas the other pulse moves toward the receiving contact with positive 
speed. The upward moving pulse reaches the injecting contact and is absorbed there before 
the downward moving pulse arrives at the other contact. In this case, the field profile close 
to the injecting contact is quasi-stationary, and a local maximum of the field is formed when 
we match this region with the trailing region of the downward moving pulse. After the 
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critical value of the current is reached, another pulse pair is nucleated, and the same process 
is periodically repeated. 

It is interesting to evaluate in some detail the process of nucleation and disappearance 
of pulses during these SSCOs. In the current vs time diagram in Fig. [3, we distinguish 
different regions depending on the number of constant solutions of Eq. (1171) that exist for 
the corresponding instantaneous value of J(t). In region A, there is only one constant 
solution on the third branch of j(F). In region B, there is one constant solution on the 
third branch and two on the second branch of j(F). In region C, there are three constant 
solutions, one on each branch of j(F), whereas two of these solutions are on the second 
branch and one of the first branch of j(F), if J(t) is in region D. There is only one constant 
solution located on the first branch of j(F), if J(t) is in region E. For a constant voltage bias, 
a pulse moving upstream may be generated only if J(t) surpasses a critical value (1.007454), 
which is located in region C and marked by an arrow. Once generated, the upstream moving 
pulses persist for any instantaneous value of the current density. These SSOC are apparently 
weakly chaotic: we have calculated the corresponding Lyapunov exponents and found that 
there is a single positive exponent with a rather small value of 7.9 xlO -6 . 

We have also observed SSCOs mediated by dipole waves that nucleate alternatively at two 
different QWs of the SL as shown in Fig. [HJ During these SSCOs, there is only one constant 
solution of Eq. (fTTj) for any instantaneous value of J(t). This solution is either on the first or 
the second branch of j(F). At about t = 3232, where J(t) reaches its global maximum, two 
dipole waves are nucleated at two different QWs. They become fully developed pulses and 
move with positive speed toward the receiving contact. When the first one arrives there, the 
current increases so that the corresponding constant solution of Eq. (TTTT) is on the second 
branch of j(F). In this case, a small dipole wave is nucleated at the local field maximum, 
where the tail of the first dipole meets the depletion layer near the injecting contact. The 
small dipole wave never grows into a fully developed field pulse, and it continues advancing, 
until the old large dipole wave disappears at the receiving contact. At this time (about 
4296), a large current spike appears, and a new dipole wave is formed closer to the injecting 
contact than the small dipole. J(t) decreases abruptly, while the small dipole disappears, 
and the newly created dipole reaches a large size and moves toward the receiving contact. 
Shortly afterward, a new current spike marks the creation of another small dipole. This 
small dipole travels toward the receiving contact, and it grows only when the only existing 
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FIG. 8: (Color online) Current vs time (top) and field distribution vs time (bottom) displaying 
dipole-mediated bulk SSCOs for an injecting contact conductivity a = 1.083425 (cf. Fig. [2]). The 
other parameter values correspond to configuration I in Table U and are N = 199, v = 8.44566, 
5 = 3.7287 x 1(T 3 , V = 1.00417485, and a 2 = 1.0523129. 

large pulse reaches the receiving contact and disappears. A small dipole formed closer to 
the injecting contact does not grow, until the large pulse reaches the receiving contact and 
disappears without triggering a new dipole wave. In this case, the corresponding pulse is 
close to the receiving contact. When it reaches the contact and disappears into it, two new 
dipoles are simultaneously triggered and become fully developed. A scenario similar to the 
one previously described follows, marked again by a large current spike. The situation is 
not repeated exactly: there are small differences in the QWs at which pulses are nucleated, 
differences in the size and lifetimes of the small dipoles, etc. These oscillations also seem to 
be weakly chaotic, in which a single Lyapunov exponent is small and positive. 

If the injecting contact conductivity is smaller so that the contact current aF intersects 
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j(F) on the NDR branch thereof, there appear standard SSCOs due to repeated dipole pulse 
nucleation at the injecting contact and motion toward the receiving contact. 

VI. CONCLUSIONS 

We have calculated the recombination coefficient as a function of the applied electric field 
for undoped, photo-excited, weakly coupled GaAs/Al^Gai-^As superlattices. Depending on 
the Al content x, the superlattice may have only one static domain for small x or two stable 
differentiated static domains for 0.45 < x < 1. In the latter dc voltage biased 

SL under weak photoexcitation may exhibit self-sustained oscillations of the current due 
to repeated nucleation of a charge monopole or dipole waves at the injecting contact and 
their motion toward the collector. For high photo-excitation intensities, the walls separating 
electric field domains are mostly pinned, and self-sustained oscillations of the current occur 
only in narrow voltage intervals. For small x among other unusual phenomena, there may 
appear weakly chaotic SSOC due to dipole dynamics in dc-voltage-biased SLs for high photo- 
excitation intensities, during which nucleated dipole waves can split in two oppositely moving 
dipoles. These dipoles are pulses of the electric field with shapes and behavior similar to 
pulses in excitable media, where a sufficiently large disturbance about the unique stable 
domain may induce them. For other parameter values, the unique static domain is unstable, 
and the underlying dynamics is oscillatory so that wave trains formed by succession of pulses 
give rise to self-sustained oscillations of the current. 
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